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ABSTRACT 

We use HSTPNFC3 imaging from the CANDELS Multicycle Treasury Survey, in conjunction with the Sloan 
Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses > 3 x 10'"Mq 
from z = 2.2 to the present epoch, a time span of 10 Gyr We explore the relationship between rest-frame 
optical color, stellar mass, star formation activity and galaxy structure. We confirm the dramatic increase from 
z = 2.2 to the present day in the number density of non-star-forming galaxies above 3 x 10'"Mq reported by 
others. We further find that the vast majority of these quiescent systems have concentrated light profiles, as 
parameterized by the Sersic index, and the population of concentrated galaxies grows similarly rapidly. We 
examine the joint distribution of star formation activity, Sersic index, stellar mass, inferred velocity dispersion, 
and stellar surface density. Quiescence correlates poorly with stellar mass at all z < 2.2. Quiescence correlates 
well with Sersic index at all redshifts. Quiescence correlates well with 'velocity dispersion' and stellar surface 
density at z > L3, and somewhat less well at lower redshifts. Yet, there is significant scatter between quiescence 
and galaxy structure: while the vast majority of quiescent galaxies have prominent bulges, many of them have 
significant disks, and a number of bulge-dominated galaxies have significant star formation. Noting the rarity 
of quiescent galaxies without prominent bulges, we argue that a prominent bulge (and perhaps, by association, 
a supermassive black hole) is an important condition for quenching star formation on galactic scales over the 
last 10 Gyr, in qualitative agreement with the AGN feedback paradigm. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — galaxies: evolution — galax- 
ies: general 



1. INTRODUCTION 

The last decade of study has brought into sharper focus 
the bimodality of the star formation histories of galaxies. 
For star-forming galaxies alone, there is a relatively tight 
distribution of star formation rates (SFRs) at a given mass 
(~ 0.3 dex scatter, with a fractio n of outlier s to high SFR; e.g., 
iBrinch mann et al. 2004, Sali m et al II2007I). pe rsisting out to 
z>2 CNoeskeet aL.2007; Wuvts et al.ll201 lal ). The red se- 
quence has SFRs substantially below those expected for star 
forming galaxies (but often with some star formation; see, 
e.g., |Yi et al. 2005). We will call these 'quiescent galaxies' 
in what follows. The relative prominence of the two popula- 
tions is a function of stellar mass, surface density, inferred 
velocity dispersion M//^ oc g^, and galaxy s tructure (e.g., 
Strateva et al.' '2001; 'Kauffmann et al.' '2003^, 'Blant on et all 
2003; Franx et al. 2008; Bell 2008; van Dokkum et al.ll20111: 



IWake et al.ll2012bt ICheung et al1l2012h . This correlation be- 
tween the structural properties of galaxies with their stellar 
populations is important: it signals that the processes that de- 
termine the structures of galaxies at least correlate, and per- 
haps are the same as, the processes that shape whether or 
not a galaxy has cold gas and star formation. Furthermore, 
these two populations evolve in their relative importance: 
while the star-forming population has a stellar mass func- 
tion that evolves slowl y (Bundv et al. 2005; Borch et al. 200^ 
iBellet al.ll2007t iPeng et al, 2010; Brammeret aL .2011), the 
quiescent galaxy stellar mass function evolves rapidly from 
z ~ 2 to the present day (largely in normalization by factors of 
~ 10, but with m odest or no evolution in shape or 'character- 
istic' mass M*; Bell et al. 2004; Borch et al.''2006; Faber et alj 
2007; Brown et al. 2007; Taylor et al. 2009; Peng et al . 201^ 
Domfnguez Sanchez eVal~i201 UlBrammer et alfcoi lb . 
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A great deal of work, both theoretical and observational, 
has been carried out to try to better understand the drivers 
of the evolution of these two populations, particularly why 
some galaxies appear capable of shutting off their star for- 
mation while others are incapable of doing so. In this study, 
we will focus on processes that can shut down star for- 
mation in galaxies which reside in the center of their own 
dark matter halo. The clear effects of gas re moval/starvation 
in dense environments (see, e.g 
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Peng et al." ^20TU, Wein mann etal] 1201 a iPeng et all 12011 
Weinmann et al. 2011 for recent discussions of this issue 
using survey datasets) are only a minor contributor to the 
evolution of 'cosmic averaged' galaxy population, owing to 
the small number of galaxies inhabiting dense environments 
dPeng et al.ll20Tot Ivan der Wei et al.ll2010i) . Accordingly, we 
do n ot discuss the eff ect of environment in this paper (see, 
e.g., Peng et al] 12010 ' for a careful discussion of the effects 
of environment as a function of cosmic time). An important 
point is that models that include the growth of the dark mat- 
ter framework, gas cooling, star formation and stellar feed- 
back alone fail to p redict a widespread population of non-star- 
forming galaxies ("Benson et al.' '2003'; i Cattaneo et al.l 120061 : 
[Somerville et al. 2008; Dave et al. 2011); all galaxies are ex- 
pected to accrete cold gas and form stars. 

A number of possible mechanisms have been proposed to 
keep a galaxy in the center of its own halo free of a signifi- 
cant cold gas content. Noting the strong tendency of quiescent 
galaxies to have prominent (or dominant) bulge components, 
it has long been thought that merging plays an important role 
in determining their struct ure (e.g.. .Tbomre & Toomre 19721 
larnes & Hernquist 1992; iHernquistil 19931; iNaab et al. 2006; 
Hoffman et ah,2010i) . There is a variety of evidence that is 



qualitatively (and in certain cases quantitatively) consistent 
with this picture: the approximate equality of the merger rate 
and the quiesce nt galaxy formation rate (e.g., Hopkins et al] 
l2010tlRobaina et al. 201 O p; the detailed kinematic structure o f 
early-type galaxies tNaab et all 120061: iHoffman et all 120101) : 
and, the empirical association between relatively younger 
stellar populations and substructure (tidal tails, shells, asym- 
metric Sj_etcJ__in_quiescen^_^^ (Schweizer & Seitzer 
[T99llTal et al 2009^. lGv6rv & Bell l'2010). 

Largely in the merger context, the possibility that feedback 
from accretion onto a supermassive black hole may drive gas 
out o f galaxies (iKauf fmann & Haehneltil2000t ISpr ingel et alj 
2005 ) or keep gas arou nd galaxies from cooling^roton et al 



2006t ISomervil le et al.ll20 08) has been explored. A wide ar- 



ray of observations are, at least at face value, qualitatively 
consistent with such a pi cture: e.g., galaxies with big bulges 
have big black holes (.Magorrian et alj [19981 : iGiiltek in et alj 
120091) . low-redshift galaxies without prominent bulges cannot 
shut off their star formation on their own (Be\] 2008), rapid 
large-scale winds are observed arou n d pos t-starburst galax- 
ies and quasars (iTremonti et alJl2007l: iProchaska & Hennawil 
120091 and the energy measured in AGN-inflated outflow 
cavities in the hot gas atmosphere of groups and clusters 
of galaxies is approxim ately consistent with the energy re - 
quired to offset coohng dBest et al.ll2006l; iFabian et al.ll2006l) . 
There are a number of other possible mechanisms for shut- 
ting off star formation, however, that are not related to feed- 
back. A few examples are: the heating of t he gaseous 
halo through virializatio n of the gas content (| Naab et al.l 
■2007; Dekel & Birnboim 120081: iKhochfar & OsFr iker 2008; 
Johansson et al.i i2009 ). changes in the mode of gas accre- 



jKeres et al.l 120051; iDeke l & BirnboimI l2006t ICattaneo et al] 
2006; Birnboim et al. 200^, and the poss ibility of the growt h 
of a stellar spheroid stabilizing a gas disk dMartig et alj2009l) . 

1.1. The goal of this paper 

Given the range of possible mechanisms for shutting off star 
formation on galactic scales for galaxies in the center of their 
own halos ('centrals' hereafter), gathering empirical insight 
into the properties of non star-forming galaxies as a function 
of cosmic epoch can be helpful. It has been argued that the key 
parameter that correlates with the paucity of star formation 
is ste llar surface density ( Kauffmann et a l. 2003; Fran x et aTl 
120081) . or possibly velocity dis persion ( Wake et al..,2012bl ; or 
roughly equivalently as M/R; Franx et al.l 120081). Ye t , for a 
sample of low-redshift galaxies from the SDSS, HI (120081) 
instead argues that Sersic index correlates much better with a 
lack of star formation for galaxies in the center of their own 
halo (that could not have been stripped of their gas by exter- 
nal influences), as non-star forming galaxies have uniformly 
high Sersic indices but a range of surfac e densities that over - 
lap with star-forming galaxies (see also ICheung et al.ll2012l) . 
Such an investigation of Sersic indices has not been carried 
out at z > 1 owing to a lack of large-scale near-IR Hub- 
ble Space Telescope (HST) imaging until recently (see, e.g., 
iKriek et al.ll2009l ICameron et aillloTol ISzomoru et aTllMlll 
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and Ivan Dokkum et al.ll201 ll for early progress towards this 
goal at z '--^ 2: IWuvts et al.ll201 Ib lexplores this in more depth). 

A not inconsiderable challenge in achieving this goal is the 
definition of what constitutes star-forming or quiescent galax- 
ies. At a given stellar mass, the SFR of galaxies depends 
strongly on redshift, evolving by a factor of 5-10 by z ^ 1 
and another factor of 4 or so out t o z ^--^ 2 dZheng et al.ll2007l : 
iNoeske et al.1 120071: iDunne et all [20091; iKarim et al.l l20Tll) . 
Noting that the la scatter in SFR at a given stellar mass 
for the vast majorit y of star-forming galaxies is ^ 0.3 dex 
dNoeske et alJ l2007l) . one could choose to define a quies- 
cent galaxy as one that has a SFR more than la below the 
star forming galaxy locus at the redshift of interest, and a 
star-forming galaxy as any galaxy forming stars at a higher 
rate. An alternative approach is to separate ga laxies by 
their optical-near-IR colors (e.g., U -V/V - J; Wuyts et al] 
120071 [Williams et al.. .2009). where galaxies dominated by 
old stellar populations are distinguishable from star-forming 
galaxies with e ven subst anti al dust reddening (as used by 
e.g.. [Wi lliams et "ani2009l and lBrammer et al]|201 ll : see also 
iPatel et al.l2 012). In this paper, we adopt both techniques. We 
note that a galaxy at z ^ 2 which is defined as quiescent ac- 
cording to these two criteria may have a SFR considerably in 
excess of almost all star-forming disk galaxies at the present 
day. While this means that our sample of quiescent galax- 
ies does not have identical properties across all redshifts, it 
does isolate a sample of galaxies with unusually low SFRs 
at that epoch given their stellar masses (confirmed by 24/im 
stacking) — one would like to understand why their SFRs are 
unusually low at that epoch. 

In this paper, we explore the evolution of the structures of 
galaxies as a function of redshift, and how they relate to the 
star formation activity in those galaxies. We use new near- 
infrared imaging from the Wide Field Camera 3 (WFC3) on 
the HST taken as part of the Cosmic Assembly Near-IR Deep 
Extragalactic Legacy Survey (CANDELS) Multi-Cycle T rea- 
sury program dGrogin et alj201 ItlKoekemoer et al.l201 ll) . fo- 
cusing on the 0.6 < z < 2.2 galaxy population in the UKIRT 
IR Deep Sky Survey (UKIDSS: .Lawrence et al. .2007.) Ul- 



Ten Gyr of quenching star formation in galaxies 



3 



tradeep Survey (UDS; These data allow us to explore 
the structure of the > 3 x 10'"Mq galaxy population 
in the rest-frame optical to z < 2.2. We supplement this 
w ith data from th e Sloan Digital Sky Survey Data Release 
2 (lAbazajian et al.l l2004) to connect with the properties of lo- 
cal galaxies. We perform two basic analyses to explore the 
evolution of the galaxy population. First, we explore the evo- 
lution of the galaxy population drawn from an 'equivalent' 
constant comoving volume as a function of redshift to get a 
sense of how the star formation and structural properties of 
the galaxy population evolve with cosmic epoch (©. Sec- 
ond, we explore the relationship between the structural pa- 
rameters of galaxies and their star formation activity using 
the full sample at each epoch (to maximize number statis- 
tics), in an attempt to understand which structural parame- 
ters best correlate with a lack of star formation activity (® . 
In what follows, we use Vega magnitudes for rest-frame col- 
ors, assume that every s tar ever fo rmed does so according to 
a universally-applicable IChabried ([2003) stellar IMF, and as- 
sume Ho = 70kms~' Mpc"', fi„,,o = 0.3 and ft\,o = 0.7. 

2. DATA 
2.1. UDS imaging data 

Our sample definition and selection is based on the pub- 
lic /T-band selecte d photometric and pho tometric redshift cat- 
alog produced by ^ Williams et al.l (I2009 f1. We adopt these 
redshifts in this paper as the basis for conversion of appar- 
ent magnitudes into rest-frame colors, magnitudes and stel- 
lar masses, and the conversion of apparent to physical sizes. 
The IWilhams et al.l (^2009) catalog use s J and .g-band data 
for the UKIDDS U DS Data Release 1 (.Lawrence et al.ll2007l; 
IWarren et al.l l200'7h in co njunction with B, R, i and z-band 
imaging from the SXDS (iFurusawa etani2008l) and 3.6^im 
and 4.5 /im IRAC imaging da ta taken as part of the SWIRE 
survey (Lonsdale et al. 2003). Total fluxes were calculated 
using an elliptical Kron ( 1980) aperture, and observed-frame 
fluxes were calculated using a matched 1"75 circular aper- 
ture on PSF-matched images (with the obvious exception of 
the poorer resolution IRAC imaging data, whose fluxes were 
measured in 3" apertures and aperture corrected to the smaller 
aperture size). Finally, Williams et al. (2009) used the HAZY 
photometric redshift code (iBrammer et al.l 12008) to estimate 
photometric redshift from the photometric catalogs. 

Spectroscopic redshifts for the UDS are relatively few in 
number and preferentially focus on brighter sources; for 
sources with Zspec ^1-1^ the redshift normalized median ab- 
solute deviation (the median absolute deviation, renormalized 
to give the same value as the RMS of a Gaussian distribu- 
tion) of Izphot - Zspec I /(I + Zspec) is ^ 0.033 with 8% catas- 
trophic o utliers ([Williams et al.ll2009h . We have also com- 
pared the lWilliams et al.l (120091) photometric redshifts against 
those of S. Wuyts et al. (in prep.), who used completely 
different (deeper) photometry and a similar photoz code to 
estimate photoz for galaxies in the CANDELS/UDS cover- 
age, finding a Az/(1 +zwiiiiams) ^ 0.055 and ~ 10% catas- 
trophic outliers (defined as having | Az|/(1 +zwiiiiams) > 0.2. 
We have confirmed that use of the photometric redshifts, k- 
corrections and stellar masses from Wuyts et al. instead of the 

' We choose to use this pubhc catalog instead of other proprietary catalogs 
to better facilitate comparison with previous works and to allow easier repro- 
duction of the results. We confirm that the results and conclusions do not 
significantly change if repeated with the redshifts, colors and stellar masses 
from the cuirently proprietary catalog of S. Wuyts et al. (in preparation). 



public Williams et al.l ( 120091) photometric redshifts plus the 
rest-frame colors and masses reported here yields no signifi- 
cant changes to our results. 

To explore the structure and morphology of the 0.6 < z < 
2.2 galaxy population, we use near-infrared F160W imaging 
from HST using WFC3. CANDELS is a HST Multi-Cycle 
Treasury program (Pis: S. Faber & H. Ferguson, FID: GO- 
12060) to image five fields on the sky using the WFC3 and 
Advanced Camera for Surveys (ACS). The CANDELS imag- 
ing of the UDS field includes 2/3 orbit in F125W and 4/3 or- 
bitsin F160W, split in to tw o epochs (see Koekemoer et ai] 
120111 and lGrogin et al.ll20TTI for more details). A total of 44 
WFC3 tiles were i maged in the UDS; w hen these are cross- 
matched with the IWilliams et al.l (12009b catalog, a total of 
0.056 square degrees is covered. 

We use also 24/im flux as a diagnostic to aid in the sepa- 
ration between galaxies with active star formation and those 
with little or no star formation. We use 24/im public data 
from the SpUDS survejflof the UDS. PSF fitting photometry 
with a 13" radius was performed, aperture corrected to total 
flux; the limiting flux of the catalog is ~ 50/iJy at the ~ 80% 
completeness (or4f7 level), and uncertainties in 24/im flux are 
of order ~ 20% (largely reflecting uncertainties in converting 
aperture to total flux, and source confusion). 

2.2. Rest-frame quantities and stellar masses 

The rest-frame magnitudes and stellar masses used in this 
paper were calculated using a set of template spectral energy 
distributions from the PEGASE st ellar population models (see 
iFioc & Rocca-Volmerangd[l997l for a description of an ear- 
lier version of this stellar pop ulation model). Such m odels 
give similar results to those of iBruzual & Chariot (2003) but 
stellar masses ^0.15 dex higher than the models by Marasto^ 
(2005£l These stellar population templates have solar metal- 
licity (as would be approximate ly expected for galax ies in the 
mass range > 3 x lO'^M©; iGaUazzi et ani2005h , and we 
fit for dust attenuation as a free parameter following ICalzettil 
(2001), with values of gas E{B-V) between -0.05 and 1.5 (to 
provide a little flexibility to fit negative attenuation to account 
for small photometry problems, etc.; in practice small nega- 
tive attenuation values are rare in the fits presented here). The 
templates include a broad range of exponentially-decreasing, 
constant or exponent ially rising star formation histories, be- 
ginning sizf (see Maraston et al. 2010 for a discussion of 
the importance of exponentially rising star formation histories 
for fitting z > 2 SEDs). 

The templates treat the evolving galaxy population self- 
consistently, in the sense that all the galaxies that are in the 
z = 2 template set also appear in, e.g., the z = 1 or z = 0.6 
template sets further along their evolutionary path. Thus, 
the galaxy population at lower redshifts is required to have 
substantial older stellar populations, driving up somewhat the 
typical age of the stars and the typical stellar mass-to-light ra- 
tio at a giv en color. This is to b e contrasted with other codes 
(e.g., Wuy ts et al.1 120081 or IPannefla et al. 2009) which have 
redshift-independent template sets, but exclude those tem- 

2 http://irsa.ipac.caltech.edu/data/SPITZER/SpUDS 

Although note that Kriek et al. (2010) demonstrate that the 
Bruzual & Chariot (2003) models appear to fit the optical-near- IR SEDs of 
galaxies with large intermediate-age stellar populations substantially better 
than Maraston (2005) models do. The overall mass scale offsets between 
j3ruzual & Chariot (2003) and Maraston (20(3) models have been discu ssed 
by a num ber of papers, and recently by e.g., IBrammer et al.l <2011l) and 
IDomfngue z Sanchez et al. (2011). 
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Fig. 1 . — Observed-frame color-color diagram for galaxies in the GOODS 
South field with spectroscopic (diamonds with black error bars) or photomet- 
lic (thin grey error bars) redshifts in the range 1.4 < z < 1.6. Overplotted in 
gray are the colors of the stars-only template (thick gray line), dust extinc- 
tion (-0.05 < E(B-V) < 1.5; included in the fitting, thin light gray lines). 
The filters were chosen to have wavelengths similar to rest-frame UVJ at the 
redshift of interest. 



plates that start star formation before the Big Bang — the typ- 
ical stellar populations of galaxies with this type of template 
fit end up being significantly younger than in the method dis- 
cussed here. It is unclear to us at this stage which method is 
more realistic; we will explore this issue with star formation 
histories drawn from a semi-anal ytic model of galaxy forma- 
tion in a future paper (see, e.g.. iLee et al.ll2010l for a version 
of this exercise at higher redshift). 

These templates were compared with photometric data 
point s of each galaxy, giv en the photometric redshift value 
from I Williams et all (120091) . and the template with the small- 
est was used to calculate rest-frame magnitudes and the 
stellar M/L ratio. Rest-frame magnitudes were calculated by 
using the SED shape of the template to work out the pre- 
dicted rest-frame magnitude of the object given the two near- 
est observed-frame bands, and then performing a weighted 
average of those two estimates of rest-frame magnitude. Stel- 
lar masses were estimated using the stellar M/Ls of the best- 
fitting template, referenced to the three longest observed 
wavelengths for overall normalization. We note that the in- 
clusion of dust, while it improves the quality of SED fit, 
leads to only modest changes in the rest-frame colors and 
stellar masses of the sample explored in this paper, given the 
overall degeneracy between the effects of dust extinction and 
stellar popu lation age in the optic al-near-IR spectral region 
(iBell & de J ong 2001). 

An illustration of this technique is given in Fig. [T] where 
we show observed-frame colors for galaxies with inferred 
stellar masses in excess of 1O'°M0 with spectroscopic (di- 
amonds with black error bars) or photometric (thin grey er- 
ror bars) redshifts in the GOODS South field (we use this 
field for this example, as it has many more spectroscopic red- 
shifts in this range; the photometric redshift sample from UDS 
shows the same trends). Overplotted with a thick gray line is 
the sequence of model composite stellar populations. Thin 



light gray lines show dust attenuation vectors corresponding 
to -0.05 < E(B-V) < 1.5. As described above, we fit these 
stellar population models with dust attenuation as a free pa- 
rameter The vast majority of galaxies are well-explained by 
the template set; note that those few that are not covered com- 
pletely by the template set are still assigned rest-frame magni- 
tudes that reflect the observed magnitudes (and therefore also 
lie off the template rest-frame colors in a similar way), as the 
observed magnitudes are used to determine rest-frame mag- 
nitudes, in conjunction with a small template-dependent cor- 
rection. 

This method, when applied to galaxies with independently- 
estimated stellar masses (from independent photometry of 
similar but not identical datase ts) in this field ( S. Wuyts et al., 
in prep.), th e GOODS-S field (iWuvts et al.l200 8). or the COS- 
MOS field (iPannella et alJl2009l) yields similar masses for in- 
tensely star-forming galaxies and masses 0.2 dex larger for 
more dusty star-forming galaxies or non-star-forming galax- 
ies (as the templates used here are more d ominated by older 
st ellar populations tha n the those used bv IWuvts et al.ll2008l 
or lPannella et al.ll2009l) . with a scatter of 0.2 dex. Rest-frame 
colors are reproduced to within 0.1 mag. When this method is 
compared with the masses of Bell et al. (2003), calculated on 
identical photometry and using a more restricted set of stellar 
population models without dust, there is no offset and 0.07 
dex scatter in stellar masses, and a scatter of less than 0.05 
mag in rest-frame U-V colors. 

We adopt the stellar mass estimates described above for the 
purposes of this paper to ensure consistency of the stellar mass 
estimation method and stellar mass scale as a function of red- 
shift (the same code was used to estimate stellar masses at 
all redshifts, and the choice of templates evolves consistently 
from redshift to redshift). We have confirmed that the system- 
atic discrepancies in stellar mass between the masses adopted 
here and those by, e.g., Wuvts et al. (2008) would operate to 
strengthen our conclusions (or in the case of the second part 
of the paper, leave them unaffected); their stellar masses for 
z > 0.6 galaxies are systematically lower by up to a factor of 2, 
with a factor of 2 scatter, and the evolution of the population 
would appear more rapid than it appears in this paper 

In what follows we adopt a mass limit of 3 X lO'^M©; the 
sample is 'complete' above this limit. Completeness, often 
meant to signify the limit above which no galaxies are miss- 
ing as a result of selection, is not straightforward to calcu- 
late for multi-band photometric redshift surveys, as the mag- 
nitude limits in the different optical/near-IR bands limit the 
recovery of photometric redshifts and stellar masses in ways 
that ar e spectral type and redshift-dependent. Willia ms et alj 
(I2009li2010 ) argue the 5a /T-band survey limit corresponds 
to a stellar mass limit of log(M5cr/MQ) 10.2; they analyze 
their sample to a limit of log(Mwiiiiam,s/A^0) ^ 10.6 to en- 
sure accurate UDS-derived galaxy size estimates (for which 
higher S/N is required). We choose to analyze the sam- 
ple to log(M/„„/M0) = 10.5. Given that completeness is 
such a challenge to calulate, we have tested this limit em- 
pirically by repeating the analyses in this paper on an in- 
dependent, currently proprietary set of magnitudes, photo- 
metric redshifts, rest-frame magnitudes and stellar masses 
drawn from deeper imaging data (S. Wuyts et al., in prep.), 
finding that our results and conclusions do not significantly 
change. The masses and rest-frame magnitudes of the sam- 
ple presented in this paper are available for download at: 
|http://www.astro.lsa.u mich.edu/~ericbell/data.php, 



Ten Gyr of quenching star formation in galaxies 



5 



2.3. Sersic profile fits of 0.6 < z < 2.2 galaxies 

To describe the structur e of th e galax y population at 0.6 < 
z < 2.2, we use parametric Sersic ( 1968) fits to the galaxy im- 
ages (A. van der Wei et al., in preparation). A surface bright- 
ness profile of the form E(r) = Tieexp[-K{j-^^" -1)] is fit using 
the GALFIT package (Peng et al. 2002), and the GALAPA- 
GOS wrapper (Barden et al, 2012). allowing the magnitude, 
axis ratio b/a, position angle, half-light radius r^, Sersic in- 
dex n, and central position to be free parameters. GALA- 
PAGOS estimates the sky value explicitly on larger scales, 
leading to more robust fits given typical image cutout sizes 
(see Haussler et al. 2007 for a detailed discussion of testing 
of our method using both simulated and deeper data; A. van 
der Wei in prep, demonstrate that the uncertainty in fit param- 
eters caused by sky estimation errors in this particular dataset 
are substantially smaller than the uncertainties that we adopt 
below from comparison of independent Sersic fits). The Ser- 
sic index n describes the shape of the light profile, where 
« = 1 corresponds to an exponential light profile and n = 4 
corresponds to a r'/"* law profile characteristic of massive, 
spheroid-dominated early-type galaxies. The Sersic parame- 
ter n is a reasonably good proxy for the ratio of bulge l uminos- 
ity to t otal luminosity, as illustrated in Fig. 14 of Simar d et alj 
(1201 ll) — systems with high n invariably host a prominent 
bulge, whereas systems with low n host a weak or no bulge 
component. At the depths typical of this imaging, uncertain- 
ties in the fit parameters are (51ogjgn O.lSdex, Sr^ ^ 18%, 
and Sb/a 0.07, as constrained from both fits of simulated 
galaxies and independent GALFITs to F125 W imaging of a 
subsa mple of z ^ 1.6 galaxies in the UDS (iPapovich et alJ 
l20T2h . 

The Sersic fits adopted in this paper are carried out on 
the F160W imaging data of CANDELS. This corresponds 
to rest-frame wavelength ranges of Arest ^ 0.55/0. 65/0. 9/im 
for z ^ 2/1.4/0.8 systems. A possible concern is that this 
change in rest-frame wavelength may affect the demograph- 
ics of the population. While the full dataset at shorter wave- 
lengths in the UDS field has not been analyzed with GALFIT, 
it is possible to test if this may be an issue using GALFIT on a 
smaller set of CANDELS F160W data in the GOODS South 
area, in comparison with published GALFITs on the F850LP 
ACS GEMS data for the extended Chandra Deep Field South 
(iHaussler et al.ll2007l) . We choose galaxies with 0.4 < z < 0.8 
for this test, where the F850LP data span the same range in 
rest-frame wavelength as the F160W data for 1.3 < z < 2.2. 
For systems with low n < I, we find a slight tendency for 
the F160W Sersic index to exceed the F850LP data (by < 0. 1 
dex), and for n>2 there is no systematic difference between 
the two sets of Sersic fits. The scatter around these modest 
offsets is ^0.2 dex, equivalent to the combined uncertainties 
of the fits. The fraction of systems with n > 2.5 in F850LP 
(restframe 0.6/im) is in fact 20% larger than the fraction 
derived using the F160W imaging (restframe ~ 1/im). This 
indicates that the fraction of n > 2.5 galaxies presented here 
at 0.6 z J; 1 -3 is likely to be close to, or perhaps up to 20% 
lower than the evolution of the « > 2.5 fraction if it were mea- 
sured in the rest frame: this operates to make the evolution of 
the population demographics more rapid still than we mea- 
sure in a fixed (red) observed band. We conclude that our use 
of F160W data alone across the 0.6 < z < 2.2 redshift range 
is not an important source of systematic error in this analysis. 



2.4. SDSS parameters for the low- redshift comparison 
sample 

In order to connect with the present-day galaxy popul ation, 
we use a sample of low-redshift galaxies e xplored in iBeH 
(2008) from the SDSS Data Release Two (lAbazaiian et al. 
I2004i) and presented in the NYU Value-Added Galaxy Cat- 
alog ([Blanton et al. 200 5). W e use foregr ound extinction- 
corrected ( Schlegel et al.lll998h . /t-corrected (Bell et al. ll2003h 
r-band absolute Petrosian magnitude for the galaxy ab- 
solute magnitude (random and systematic uncertainties < 
0.15 mag), and model colors for higher S/N es timates of 
galaxy color (uncertainties < 0.05 mag). Following lBell et al.l 
(l200l . we have merged this catalog with the Two Micron 
All-Sky Survey (2MASS; Ski-utskie et al. 2006) to facilitate 
SED fitting and to allow splitting of galaxies by U -VIV -J 
into quiescent and star-forming populations. Stellar masses 
and rest-frame colors were estimated from ugrizJK photome- 
try using the above stellar population model templates. 

Star formation and AGN classification s and estimates of 
total star formation rate were taken from 'Brinchman n et alJ 
(12004) using emission line measurements described in 
iTremontfet al. (2004). Galaxies are classified as star forming, 
AGN, composites, or are left unclassified (typically because 
the galaxies lack line emission in their SDSS spectra). 

The only source of Sersic fits for all galaxies in our sam- 
ple is Blanton et al. (2003fl Blanton et al. (2003) fit the light 
profile of galaxies in the SDSS, measured in circular aper- 
tures, with a seeing-convolved Sersic (.1968i) pr ofile for all of 
the ga laxies in the VAGC. The Sersic fits by iBlanton et alj 
(12003b give values for the Sersic index, in particular, that 
are offset from n values determined using 2D galaxy im- 
age fits. We have comp ared the Sersic indices (and other 
fit parameters) f rom B l anton et al.l (2003) to fits carried out 
by Ivan der Well (l2008h on a small subset of galaxies in the 
NYU VAGC. We find that the Sersic index estimates are 
related: logion2D -0.39 + 1.75 logiQnNYUvAGC, with 0.2 
dex scatter Half-light radii show the following correlation: 
logio'"2D ^ log,„rNYUVAGC-0.05 + 0.025n2D,new, with 0.1 dex 
scatter, where niD.new is the estimate of equivalent 2D Ser- 
sic index derived from the NYU VAGC Sersic index using 
the above relation. A similar analysi s was carr i ed out with 
completely independent GALFITs by ' Guo et a D (l2009h . and 
importantly the above trends are identical in the case of «, 
and < 0.05 dex different in the case of r, to the median offsets 
as a function of n in their Fig. Al. Recall that the 0.2 dex 
scatter between th e ' rescaled' NYU n values and those of 
Ivan der Well (120081) or lGuo et al.1 (120091) (or the 0. 1 dex scatter 
in radii) is comparable to the typical joint uncertainties in any 
comparison of even unbiased values of n (or r). We conclude 
that these rescaled NYU n and r values are unbiased, have un- 
certainties comparable to those determined directly from 2D 
fits, and are appropriate for connecting the evolution of n and 
r with the results of 2D fitting for the z > 0.6 galaxies. 



^ Si mard et alj (201 Ih also fit all galaxies in the SDSS with single Sersic 
profile fits. Owing to the low redshift of our SDSS sample, about 1/8 of the 
galaxies in this particular sample lack fits in Simard et al. (201 1), as they are 
above the bright limit adopted for the construction of those catalogs. Accord- 
ingly, we do not adopt the estimates of Simard et al. ( 20 11 ) for this work. The 
vast majority of galaxies that do have fits agree well with the trends reported 
below between NYU and Sersic fits performed directly to the imaging data 
(put differently, connected NYU fits agree with Simard et al. 2011). Results 
would be unchanged if we adopte d thos e estimates instead for the fraction of 
galaxies with fits in Simard et al] 120 llh . 
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3. THE EVOLVING RELATIONSHIPS BETWEEN STELLAR MASS, 
COLOR AND MORPHOLOGY 

3.1. Separating galaxies into quiescent and star fanning 
populations 

We separate galaxies by their star formation activity using 
two sets of independent diagnostics: mid infrared (or, for the 
SDSS, emission line) information, and position on optical- 
near-infrared color-color diagrams. We use both cuts in this 
paper; we describe the 24/im-derived cuts first (showing how 
they relate to rest-frame color cuts) and then show the rest- 
frame color cuts (showing how they relate to 24/im cuts). 

At z > 0.6, galaxies are classed in large part according to 
their 24/im emission properties as an admittedly imperfect 
proxy for obscured SFR (we use emission-line diagnostics 
and SFR estimates from Brin chmann et al.ll2004l for galaxies 
from the SDSS). Substantial 24/im flux can also result from 
AGN activity. We do not attempt to discriminate between 
AGN activity and SF activity for our purposes here, simply 
noting that systems with 24/im flux dominated by AGN at this 
24/im lumin osity and redshift range are not the dominan t pop- 
ulation (e.g.. iDonlev et al.l 120081: iKartaltepe et al.ll2010 ) and 
that we are primarily attempting to weed out galaxies whose 
rest-frame optical colors are a poor reflection of the stellar 
populations in that galaxy, a goal for which our simple ap- 
proach is sufficient. 

In this spirit, we wish to avoid an e xplicit, and uncertain , 
conversion of 24/im flux into SFR (e. g., lPapovich et al.ll2007l : 
lElbaz et al.ll20Tlt IWuvts et alj|201 lah . At each redshift of in- 
terest, we fit the relationship between 24/im flux and stellar 
mass for galaxies detected at 24/im (shown as filled symbols), 
as shown in Fig. [D The approximate trend at all redshifts is 
logm/24//iJy ^ 2 + 0.51ogm(M*/3 x lO'^Mg), with a scat- 
ter of less than 0.3 dex (we use the actual fits and scatters, 
which vary slightly with redshift, to perform the split into star 
forming and quiescent). The slope and scatter of this relation- 
ship is well-documented and studied (e.g., Salim et al. 2005; 
IZheng et alJl2007t iNoeskeet^ilfMOTt iKarim et alJl201 Ih : it 
is a remarkable coincidence that the zero point in terms of 
24/im flux varies so little with redshift, owing to the inter- 
play between the dramatic reduction of SFR at a given stel- 
lar mass with decreasing redshift, the luminosity distance and 
the redshift-dependent 24/im A:-correction. Quiescent galax- 
ies then must have a 24/im flux (UDS) or SFR (SDSS) lower 
than -1(7 from the star forming galaxy locus. 

In Fig. 121 we have color-coded symbols by their position on 
the rest- frame U-V/V-J diagram (Fig .[3]l, using the slightl y 
redshift-dependent cuts described in I Williams et alj (120091) . 
Orange symbols show galaxies with rest-frame optical-near- 
IR colors characteristic of quiescent galaxies, and purple sym- 
bols show galaxies with colors characteristic of star-forming 
galaxies with a range of reddening values. In Fig. [3] we show 
the optical-near-IR colors of galaxies in the six redshift in- 
tervals of interest. In contrast to Fig.|2] we have color-coded 
the symbols in Fig.[3]by 24/im flux (UDS) or SFR (SDSS). 
Galaxies with 24/im fluxes/SFRs lower than -1 cr from the star 
forming galaxy locus have been color-coded red, and galaxies 
with fluxes/SFRs higher than -la from the SF galaxy locus 
have been color-coded blue. In the z < 0.05 slice, emission- 
line diagnostics are available, and any object with AGN-like 
lines or composite SF/AGN lines (Brinchmann et al. 2004) 
has been color-coded green. 

Inspection of Figs. |2] and [3] shows the large degree 
of overlap, and the complementarity of having both ex- 



plici t 24/im/SFR inform a tion and U — V/V — J color s (see 
also IWilliams et al.1 120091 IWuvts et all 12009. Bramm er et al.l 
and a morphological investigation by P atel et al.ll20I2h . 
Galaxies with quiescent U - V/V - J tend, for the most part, to 
be undetected at 24/im (Fig.|2|i. There are exceptions to this: a 
few galaxies with quiescent U-V/V-J at z < 1 .5 are detected 
at 24/im, and by z 2 it is clear that for our particular dataset 
the contamination of the quiescent region of U - V/V - J 
color space by 24/im-detected objects is significant. Stack- 
ing at 24/im of the remaining individually-undetected galax- 
ies with quiescent U-V/V-J yields marginal {2-3a signifi- 
cance) detections at the 5- 10/iJy level at all redshifts (indi- 
cating SF/AGN activity a factor of > 10 lo wer than typical 
star-f orming galaxies at that redshift; see also iPapovich et all 
'2006^ for an early discussion of star formation in red-selected 
galaxies; these measurements are not shown as they fall off of 
the range of data values plotted). 

Conversely, galaxies detected clearly at 24/im are almost 
always in the star forming region of U-V/V-J (Fig. [3]), 
but again with some exceptions (e.g., at z 1.4 there is a 
clear group of galaxies with star forming colors that are in- 
dividually undetected at 24pm). Stacks of those few indi- 
vidually 24/im undetected galaxies with star-forming colors 
yields significant detections at the 15-25/iJy level (blue filled 
points with error bars on Fig.|2]i, a factor of a few below the 
SF galaxies locus, largely consistent with an interpretation of 
these systems as the low SFR tail of the SF galaxy population. 
This high degree of correspondence between the two meth- 
ods has been shown before by e.g.. lWilhams et al.l (2009) and 
IWuyts et"al] (2009). Fig. [3] also shows that the use of rest- 
frame color information for z < 0.05 is particularly valuable; 
galaxies classified as AGN can have either quiescent or star 
forming colors. 

We separate galaxies into quiescent and star forming us- 
ing both criteria to capitalize on their different strengths and 
shortcomings. Discrimination by U-V/V-J is sensitive to 
lower amounts of star formation than lAfjm separation, es- 
pecially at z ~ 1.4 and z ^ 2. On the other hand, sepa- 
ration by 24/im is considerably less sensitive to photoz er- 
ror than U - V/V - J, as one simply needs to know which 
redshift bin the galaxy is in, and even some cross-talk be- 
tween bins can be tolerated. We define quiescent galaxies 
as having both 'quiescent' colors in U-V/V-J and 24/im 
fluxes/SFRs lower than - la from the star forming galaxy lo- 
cus at the redshift of interest. For SDSS galaxies with AGN- 
like emission lines, we split only on the basis of U-V/V-J . 
Star forming galaxies are defined as those which satisfy ei- 
ther (or both) of the U - V/V - J star forming galaxy color 
cuts or having 24/im fluxes/SFRs brighter than - la from the 
star forming galaxy locus. One can see that incorporating 
24/im data and insight from U-V/V-J into this analysis is 
crucial. As can be seen directly in Figs. [2] and [3] and fur- 
ther appreciated by the intermixing of blue and red symbols 
at red rest-frame tZ-V in Figs. |4H8] failure to flag galaxies 
by multiwavelength-derived SF activity leads to considerable 
confusion between star forming galaxies with substantial dust 
columns and non star-formi ng galaxies, esp ecially at higher 
redshifts (ITaylor et al.ll2009t iBrammer et al.ll20I II) . 

3.2. Evolution of the galaxy population in a given comoving 

volume 

Fig. |4] shows the rest-frame U -V colors of galaxies, as a 
function of their stellar mass, in six different redshift bins. We 
choose to show the properties of the galaxy population as a 
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Fig. 2. — Flux at 24fim as a function of stellar mass in three different redshift bins; detections are shown as filled symbols and 4(t upper limits shown by 
arrows. The thick solid line and dotted lines are the robust linear fits to the data and the ilcr range. Colors denote classification according to U — V/V -J (Fig. 
[5): purple symbols denote galaxies that lie in the star forming part of the color-color diagram, and orange galaxies lie in the quiescent region of the plot. The 
large solid purple points with error bars show the stacked 24/im flux and uncertainty for galaxies individually undetected at 24fim but in the star-forming part of 
the U - V/V -J diagram. 
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Fig. 3. — Rest-frame U — V as a function of V-/ color for the six redshift 
bins used in this paper. Galaxies are color-coded according to their 24/^m 
properties (z > 0.6) or emission line properties {z. < 0.05); (quiescent) galax- 
ies not detected or fainter than -Icr from the star forming galaxies locus are 
color-coded red, and the rest of the galaxies (all star forming) are color coded 
blue. Green galaxies at z < 0.05 are galaxies in the SDSS with AGN-like 
emission lines. Superimposed are the rest-frame color cuts used in this paper, 
following Williams et aL 12009) . All galaxies have stellar masses in excess 
of 3 X IO'^Mq. 



function ofU-V re st-frame color to connect with other stud- 
ies (e.g., iBell et all [2004; .Borch et al.. .2006: .Ruhland et al.. 



I2009t IWhitaker et all IMa iBrammer et alJ 1201 ll) and as a 
joint (rather sensitive) constraint on SFH and dust content. 
Galaxies are color-coded by star formation activity ( ^3.11 ): 
red symbols show galaxies classified as quiescent using both 
24/im information and U-V/V-J colors, and blue symbols 
show the remaining population. In all panels of this figure, the 
galaxy population has been Monte Carlo subsampled down 
to an equivalent comoving volume of lO^Mpc^ by adjusting 
the number of galaxies to track the number density of galax- 
ies with > 3 X 1O"^M0 determined from larger surveys 
(the line in Fig. \T3[ . Put differently, variations in the num- 
ber of galaxies from panel to panel illustrate true evolution in 
the galaxy population (as the volume is fixed; see Appendix 
|A]for further discussion). Filled symbols show galaxies with 
n > 2.5. Open symbols show galaxies with « < 2.5. In all pan- 
els, the linear size of the symbol scales with (1 -l-logjo re/kpc), 
where is the half-light semi-major axis, and the axis ratio 
of the symbol is the same as that of the galaxy of interest. The 
black line is shown in all panels for reference at the approxi- 
mate position of z ^ 2 non star-forming galaxies. 

The evolution of the galaxy population in the epoch z ^ 2 
to the present day is obvious. As has been argued by a nun> 
her o f other author s (e.g., Arnouts et al. 2007; Fontana et 
2009'; 'Tavlor et al.' '2009; Ilbertetal. 2010; Cassata et 
2011: DomJnguez Sanchez et al. 201 1 ; Brammer et al. 201 1| ) 



there is dramatic evolution in the overall number of galaxies 
with > 3 X 1O'°M0 (as quantified in Fig. [131). Further- 
more, Fig.|4]shows that the evolution of the number density of 
quiescent galaxies is particularly striking (again, as has been 
argued by the above cited works). Fig. |4]makes it clear how- 
ever that the evolution of the star formation activity of the 
intermediate-mass and massive galaxy population is accom- 
panied b y a large-scale ch ange in the structure of galaxies 
(see also'Wuvts et al.ll201 l b, who see similar behavior) from 
a z > 1.5 population dominated by low n (little or no bulge), 
mostly star-forming systems to the present population, domi- 
nated by galaxies with high n (with a prominent bulge), many 
of them quiescent (but not all of them). Quantitatively, there 
are a factor of ~2.5 (3.5) more quiescent (n > 2.5) galaxies 
with M^, > 3 X IO'^Mq today than there were galaxies with 
those masses at z ~ 2, respectively. Thinking about it differ- 
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Fig. 4. — The evolution ofU-V rest-frame color (in Vega magnitudes) as a function of stellar mass, in six different redshift bins. Galaxies that 
appear in this figure are subsampled to a fixed comoving volume of lO^Mpc' at all redshifts to show the evolution of the massive galaxy content 
of a 'representative' volume with cosmic time. Open symbols show galaxies with n < 2.5 and filled symbols show galaxies with n > 2.5. For 
all symbols the axis ratio of the galaxy is reflected by the axis ratio of the symbol, and the size of the symbol scales with (1 +log,(, r^/kpc). The 
black line is shown in all panels at the approximate locus of z ~ 2 red galaxies. Galaxies are color-coded by star formation activity — galaxies 
classified as quiescent by both 24/im/SFR and U -V/V -J diagnostics are color-coded red, and blue symbols show all other galaxies. AGN 
with z < 0.05 are classified using only U - V/V- J. 
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ently, the current population of quiescent (« > 2.5) galaxies 
with M, > 3 X 1O'*'M0 is approximately as numerous as the 
entire M, > 3 x IO'^Mq population at z ^ 1.1 (0.7) respec- 
tively. This change in global demographics from z ^ 2 to the 
present day makes it clear that, in addition to processes that 
shut off star formation on galactic scales, there must also be 
(the same or different) processes that lead to an associated 
change in the surface brightness profiles of galaxies over the 
same time period (and given the correspondence between a 
lack of star formation and structure, the time scales of such 
processes must be comparable). 

Fig. |4] shows also that the scatter in the quiescent 
galaxy color-magnitude relation (CMR) decreases towards 
lower redshift. The evolution of CMR scatter from z ^ 
2 to the present day is well-documented in the literature 
dRuhla nd et al. 2009; Whitaker et al. 2010). The scatter in the 
U-V colors of non star-forming galaxies with U -V > 0.6 at 
z > 1.3 is 0.17 mag (our measurement) , very consisten t with 
the carefully-measured results of Whit aker et al.l (120101). who 
find a scatter of 0.13-0.2 mag for 1.3 < z < 2. iRuhland et al.l 
(|2b09) find that the scatter inU -V color at z < 1 is 0. 1 mag 
(measured much more carefully than the CMRs presented in 
Fig-El our measurements also give a scatter of 0. 1 mag), es- 
sentially indep endent of redshift. M odeling presente d in both 
IRuhland et all (l2009l for z < 1) and lWhitaker"etan (l20Toi at 
1 < z < 2) shows that the evolution of CMR scatter is natu- 
rally interpreted as being caused by a constant inflow of new 
galaxies onto the red sequence at the observed number density 
growth rate. 

4. EMPIRICAL CORRELATIONS BETWEEN A LACK OF STAR 
FORMATION AND GALAXY STRUCTURE 



4. 1 . Broad trends 



One of the most no tab le tren ds s een by lArnouts et al.l 
J2b07), iTavlor et all (l2009h Jllbert et aD J20TC 
iDommguez Sanchez et al.l (1201 ll) . iBrammer et al. (2011) 
and elaborated upon in Fig. |4] is the dramatic growth of 
the quiescent galaxy population from z ~ 2 to the present 
day. Fig. |4] demonstrates also that with the growth of the 
quiescent galaxy population comes a concurrent growth of 
the population of c o ncentrated n > 2.5 galaxies (see also 
IWuyts et al.1 HoHB- iFranx et al.1 (l2008h argues that these 
galaxies also have high surface density and M /r^, which 
should scale with velocity dispersion, and that M /r^ is the 
parameter that best correlates with a lack o f SF activity at 
z ^ 2 (see also ! Wake et al.ll2012bl:ICheun g et al. 2012). In this 
section, we explore how the different structural parameters 
correlate with star formation activity in an attempt to gain 
possible insight into the processes that drive galaxies into 
quiescence. 

In Fig. |5] shows the color-mass trends using the full UDS 
and SDSS DR2 datasets, as opposed to sub-sampling down 
to an equivalent volume of lO^Mpc^, in order to delineate the 
trends with better fidelity than the subsample shown in Fig. ID 

In Fig. |6l we show rest-frame color varies with veloc- 
ity dispersion estimated from stellar mass, half-light radius 



and Sersic index, where: a = 



CM, 



and K^{n) = 



0.954+ ■ 



73.32 



K^(n)re^(hlar 

10 465+(„_o 94)2 scales M/ re in a physically-motivated 
way t o account for the s tructure of a galaxy via the Sersic 
index dBertin et alj|2002l amounting to velocity dispersions 
at a fixed M/r^ that are ~ 0.1 dex lower for n <2 systems 
compared to those with n ^ 4). Such a scaling permits re- 
covery of observed velocity dispersions of galaxies in the 



SDSS as a function of ph otometric paramete r s to an accuracy 
of~ 0.12 dex (see, e.g.. iTavlor et al.l I2OIOI iBezanson et alj 
I2OI Ih . This velocity dispersion estimate scales also with rota- 
tion velocity for rotationally-dominated systems, albeit with a 
different proportionality constant. 

In Fig.|7] we show the variation in rest-frame color with sur- 
face de nsity within the half-hg ht r adius S = O-pM/n r^ (fo\- 
lowing iKauffmann et al.l 12003' and 'Fra nx et al] 1200 8', and as 
explored for low redshift by Bell 2008). In all figures, sym- 
bols are coded by Sersic index, axis ratio, size and star for- 
mation activity as in Fig. |4] with the exception of the SDSS 
sample where the number of galaxies allows color-coding by 
star formation activity alone. 

In Figs. [4l— H one can see tre nds previo u sly re - 
ported by Kauffmann et alj ( l2003h . iFranx et alj (120081) . 
van Dokkum et al. (201 1) or many subsequent studies: galax- 
ies with high stellar mass, high velocity dispersion, or high 
surface density, tend not to form stars (where the latter study 
is particularly relevant owing to its use of HST-deiived struc- 
tural parameters and the equivalent width in Ha as a star for- 
mation indicator). Yet, one can see also evidence that M*, 
cr' or S fail to give a complete picture of which galaxi es are 
quiesce nt (see also Bell 2008, Cheung et al. 2012, Wake et aP 
HoilH). A significant fraction of low stellar mass galaxies are 
quiescent (therefore stellar mass is a relatively poor predictor 
of quenching), and a small fraction of galaxies with z < 1 .5 
and intermediate or low values of a' and S are quiescenfl 

Fig-Eshows the trend in rest-frame color with Sersic index 
(recall that Sersic index corr elates with the rela tive promi- 
nence of a bulge component: ISim ard et al. 1201 lb . At z > 1, 
this trend was not explored by Fran x et alT i 2008h . as they did 
not analyze large-scale HST imaging, and therefore lacked re- 
liable measurements of surface brightness profile shape; this 
can, however, be compared with (and is consistent with) Fig. 
1 of lWuytset"aD(l2011b) . Symbols are similar to previous 
figures except that filled/open symbols now denote galaxies 
with above/below the median a' at that redshift for the galax- 
ies with > 1O'"'^M0 in our sample (recall that the previ- 
ous filled/open distinction was by Sersic index, which would 
be redundant). Focusing on the z > 0.6 points, one can see 
that galaxies with high Sersic index are much more likely to 
be non star-forming than their low Sersic index counterparts. 
Furthermore, one can see that there is a range of velocity dis- 
persions at a given Sersic index. A similar qualitative behav- 
ior is seen for the SDSS galaxies, notwithstanding quantita- 
tive differences in the definition of Sersic index for the NYU 
VAGC catalog. 

4.2. Which parameter correlates best with a lack of star 
formation? 

Figs.|3-l8]demonstrate that 'typical' quiescent galaxies have 
higher mass, 'velocity dispersion', surface density and Sersic 
index than 'typical' star-forming galaxies. In this section, we 
explore further which parameter correlates the best with a lack 
of star formation activity. 

Fig. |9] shows the quiescent fraction (evaluated in running 
bins of 101 galaxies at z > 0.6, or 501 galaxies at z < 0.05) as 

The actual fraction of galaxies forming stars at low E or cr' may be rather 
higher, as the sample is hmited by stellar mass. Fig.|4]shows that quiescent 
galaxies are smaller at a given stellar mass than star-forming galaxies, there- 
fore it is possible that if Figs. |6] and [7] were a' or E limited samples they 
would show a more prominent population of (lower stellar mass) quiescent 
galaxies with relatively low cr' or S. 
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a function of the rank of a galaxy in stellar mass (dotted line), 
'velocity dispersion' (grey solid line), surface density (dashed 
line) and Sersic index (black solid line) in three broad redshift 
bins. At all redshifts, stellar mass is a poor predictor of qui- 
escence. At z < 0.05, Sersic index is clearly a better predictor 
of quiescence than any other parameter; in particular, galaxies 
with low n overwhelmingly host detectable star formation. 

We quantify this by introducing the quantity Aq, 2-0.6, which 
quantifies the fractional difference in rank between the galaxy 
population being 20% quiescent and 60% quiescent (i.e., the 
rank difference corresponding to when the lines cross quies- 
cent fractions of 0.2 and 0.6). At z < 0.05, Ao.2-0.6 is un- 
defined for stellar mass, surface density and velocity disper- 
sion (as the quiescent fractions never go below 0.2), and is 
0.31 ±0.03 for galaxies when ordered by Sersic index. Put 
differently, when one orders the galaxies by Sersic index from 
low to high, from the point when the quiescent fraction is 0.2, 
one needs to go through 3 1 % of the galaxies to reach the point 
where the quiescent fraction reaches 0.6, and that bootstrap- 
ping of the galaxies being used to calculate this quantity leads 
to a ±3% variation of Aq 2-0.6- 

At 0.6 < z < 2.2, one can see a rather different situation: 
our estimate of velocity dispersion, surface density, and Ser- 
sic index all correlate comparably well w ith SF activity (the 
first and last trends agree with Fig. 4 of Ivan Dokkum et aLl 
l20Tlh . For 0.6 < z < 1.3, Ao.2-0.6 is 0.31 ±0.06, 0.40±0.09 
and 0.47 ±0.07 for n, a' and E respectively. For 1 .3 < z < 2.2, 
Ao.2-0.6 is 0.30±0.07, 0.30±0.06 and 0.37±0.06 for «, a' 
and E respectively. Echoing in a muted way the behavior 
of the SDSS sample, one can see that the Sersic index still 
corr elates with quiescence w ell (and as well as z < 0.05; see 
also 'van Dokkum etalll2011h . The correlation of star forma- 
tion activity with a' is as strong as that of Sersic index at 
1 .3 < z < 2.2, and appears to weaken with decreasing redshift. 
The correlation of star formation activity with E is marginally 
poorer than with Sersic index at z > 0.6. 

The small numbers on the lower left and right hand cor- 
ners of each panel of Figs. |5]-[8] also help to illustrate this 
point. These numbers show the fraction of quiescent galax- 
ies for two different subsamples: the half with lowest stellar 
M,, a', E, or «, and the half with highest M,, a', E, or «. We 
note that the quiescent fractions split by mass are qua ntita- 
tively similar to those in Fig. 4 of lBrammer et al.l(l201 Ih . The 
Sersic index is the metric that tends maximizes the contrast 
between the two halves of the sample (except at 0.9 < z < 1.3 
and 1.5 < z < 1.8, where Sersic index still correlates very 
well with quiescence). Repeating these analysis with alternate 
photometric redshifts, stellar masses and rest-frame properties 
(S. Wuyts et al., in prep.) produces very minor changes, with 
a slightly stronger preference for n as the parameter that best 
correlates with quiescence. 

Given that n and a' correlate well with quiescence, it is in- 
teresting to explore star formation activity as a function of 
both parameters. In the same three broad redshift bins, we 
show the U -V color of galaxies as a function of a' . Sym- 
bols are coded as they are in Figs. [4^7] The top panels show 
galaxies with n < 1.5 (galaxies with little or no bulge) and 
the bottom panels show galaxies with « > 3 (galaxies with a 
prominent bulge). The population with n < 1 .5 is overwhelm- 
ingly star forming. In strong contrast, the n > 3 population has 
a large quiescent fraction. There is a strong tendency for high 
« galaxies to have high a' (although there is considerable scat- 
ter in velocity dispersion at a given Sersic index), and not all 
high velocity dispersion or Sersic galaxies lack star formation. 



Neither parameter perfectly predicts quiescence, although it is 
clear that both perform very well at z > 1, and Sersic predicts 
quiescence better for local samples. 

There are two emergent themes that we wish to draw the 
reader's attention to. First, Figs.[8]and[T0]show that, with very 
few exceptions, galaxies with low Sersic index all appear to 
form stars at all z < 2.2. The threshold appears to be some- 
where around n ^ 1.5-2: at n < 2, the fraction of quiescent 
galaxies is < 10% (and in many redshift bins it is less than 
a few percent). When investigated in more detail using the 
same SDSS sample (Bell 2008), it was found that i) real low 
n quiescent galaxies are all satellite galaxies in galaxy clus- 
ters, i.e., they are stripped disk galaxies, and ii) the few qui- 
escent 'low «' systems in the centers of their own halos that 
remained were in fact the result of measurement error in n, 
as visual inspection showed a distinct bulge component. We 
show examples of some n < 2 quiescent galaxies at z > 1.5 
in Fig. [TT](at 0.9 < z < 1.5 there are only three « < 2 qui- 
escent systems, and those look similar to the z > 1 .5 exam- 
ples; the 0.6 < z < 0.9 points are bad Sersic fits). While some 
systems are relatively extended and have low «, and the one 
inclined galaxy is clearly reminiscent of a disk, most appear 
spheroidal and compact. Given Figs.[8]and[TT]in concert, it is 
clear the vast majority of quiescent g alaxies have a prominent 
spheroid. This extends the results of iBelll (120081) determined 
for nearby galaxies and lCheung et al.l (120121) for z ~ 0.65 to 
z ^ 2.2, when the Universe was ^ 1/4 of its present age: 
galaxies lacking a prominent bulge appear to have great diffi- 
culty shutting off their own star formation on galactic scales. 

Second, Figs. Isl and [Tol make it clear that having a high 
Sersic index alone (or indeed, having high n, a', E and M,) is 
not enough to ensure a lack of star formation. At all redshifts, 
a small minority of high n sources form stars at an appreciable 
rate. This illustrates a key point of this paper; it appears that 
for all z < 2.2 a large bulge is necessary to stop star formation, 
but is not suff icient to sto p star formation. This extends the 
conclusion of iBelll ( 120081) determined for local galaxies out 
to z 2.2, when the non star-forming galaxy population was 
considerably less prominent. 

5. DISCUSSION 

There are two main observational results in this paper: the 
rapid growth of the quiescent galaxy population between z = 
2 and the present day, and the recognition that this growth 
appears to be intimately linked to the growth of galaxies with 
prominent bulges (as quantified by high Sersic index, inferred 
velocity dispersion and surface density). 

5.1. Musings on the mechanisms that prevent significant cold 
gas in galaxies 

These results have some bearing on understanding which 
mechanisms lead to quenching of star formation in galaxies. 
Recall that the role of environmental quenching is relatively 
minor in our 'cosmic averaged' population evolution, and that 
we are focusing on which types of physical process lead to 
quenching of star formation i n galaxies in the c enters of their 
halos (the mass quenching of iPeng et al.ll20ldl) . To facilitate 
this, we will set up two straw person hypotheses: suppression 
of star formation by feedback (either star form ation or AGN 
feedback; 'feedback quenching'. e.g.JKauffma nn & Haehnel3 
I2OOOI ICroton et al.ll2006i fSomerviUe et al.ll2008,) . or suppres- 
sion of star formation because the halo reaches a certain crit- 
ical mass ('halo quenching'; e.g.. iDekel & BirnboimI l2006l 
,Cattaneo et al...2006.) . 
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Fig. |9] demonstrates that internal properties (Sersic index, 
and at z > 1 also 'velocity dispersion' and surface density) 
correlate strongly with quiescence, whereas galaxy mass only 
weakly correlates with star formation activity . Galaxy mass is 
expected to correlate well with halo mass. More et al.l (|2009|) 
measure the s catter in lum inosity at a given halo mass to be 
0.16dex, and I Yang etaH (2009) measure the scatter in stel- 
lar mass at a given halo mass to be 0.17dex, both substan- 
tially less than the ~ Idex dynamic range probed in this work. 



The observed weakness of correlation between quiescent frac- 
tion with stellar mass, coupled with the expected modest scat- 
ter between stellar mass (via its proxy, luminosity) and halo 
mass, implies a weak correlation between halo mass and qui- 
escence at all z < 2.2. This is consistent with the claim by 
iMore et akl (1201 Ih that at fixed stellar mass there is no differ- 
ence between the average halo masses of quiescent and star- 
forming central galaxies. On this basis, it would be natural 
to conclude that quiescence is not determined by halo mass 
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alone. 

We caution that such a conclusion may be premature. 
I Wake et al](l2012al) find that at high velocity dispersion/stellar 
masses the clustering of galaxies (a reflection of character- 
istic halo mass) is better described as a function of veloc- 
ity dispersion than stellar mass. On this basis, our finding 
that quiescent fraction is quite strongly correlated with veloc- 
ity dispersion may suggest that quenching is a strong func- 
tion of halo mass. We not e, however, that the sample stud- 
ied by iWake et aTl ([2012a') is complete only at high mass, 
> 1O"M0, where iMore et al. (2009) infer mcreasmg scat- 
ter in h alo mass as luminosity increases. T aking Mor e et al] 
(l2009h . lYang et al.1 (|2009.) and IWake et al.1 ([20 12a) together, 
it appears possible that luminosity correlates well with halo 
mass at lower luminosity/stellar mass and that velocity dis- 
persion correlates better with halo mass at higher luminos- 
ity/stellar mass. This issue cannot be resolved here, and the 
issue of how the weakness of the quiescence-stellar mass cor- 
relation should be interpreted by necessity remains open. 

The results here demonstrate that systems with high Ser- 
sic index, 'velocity dispersion' and high sur face density are 
much more likely to be quiescent (see also Franx et al. 2008, 
Ivan Dokkum et al. 2011, Wuyts et al. 2011b, Wake et al. 
|2012b). All of these metrics correlate strongly with the 
relative prominence of a bulge component. Given correla- 
tions between black hole mass, bulge mass, velocity dis- 
persion an d Sersic i ndex ("Haring & Rixl l2004t iPeng et al.l 
|S)06; Gra ham & Dri ver 2007; GultekiiTetdJIZOQgl), our^ 
suits tentatively support the AGN feedback paradigm, at least 
at the qualitative level. In particular, the seeming inability 
of galaxies with low Sersic indices/inferred velocity disper- 
sion/surface density to shut off their star formation is very nat- 
urally interpreted in this framework (Bell 200^ — no super- 
massive black hole, no shut-off of star formation. We note that 
feedback from star formation-driven winds may also be of rel- 
evance, but we caution that the results shown here argue that 
such winds must be efficient at wholesale ISM removal only 
during star formation events that create a prominent bulge (as 
disk galaxies keep forming stars with enthusiasm). 



5.2. Quiescent disks? 

There has been much recent discussion of the pres- 
ence and importance of stellar disks in quiescent galaxies. 
At z > 1.5, Stockton et al. ( 2004), McGrath et al. (2 008,), 
Ivan Dokkum et al.1 (l2008F and Ivan der Wei et al.1 ( 1201 lb have 
argued that most quiescent galaxies with masses in excess 
of 1O''M0 have prominent stellar disks. iBundv et aTl (120101) 
have explored this issue at z < 1 .2, finding that a large frac- 
tion of quiescent galaxies have disks in addition to signif- 
icant bulge components; indeed, the existence of SOs and 
disky ellipticals in al l environments is well-known (see, e.g., 
Ivan den Berghl l2009l and references therein). The highly- 
inclined fraction of these systems are clearly visible in Fig. 
[8] as the elongated symbols with low b/a. Such systems ap- 
pear to be somewhat more common at z > 1 .5, but are present 
(especially at lower Sersic indices n ~ 2-3) at all redshifts. 

There are two comments that we wish to make about qui- 
escent disks. Firstly, the vast majority of th ese systems 
have Sersic indic e s n > 2. In iMcGrath etaTI (2008) and 
Ivan der Wei et al.l (1201 II) . bulge/disk decompositions were 
carried out, and showed that these systems with relatively high 
Sersic index are also well-explained as composite bulge/disk 
systems with relatively large bulges. Furthermore, kinematic 
studies of local quiescent galaxies have demonstrated that the 
vast majority of quiescent galaxies have significant rotation 
(jEmsellemet al. 2011), and the incidence of strong rotation 
signatures in quiescent galaxie s does not change from z ~ 1 
to the present day (at ^ 60%, Ivan der Wei & van der Marell 
120081) . The picture that emerges is that the vast majority 
of quiescent systems have undergone some event that both 
steepens their light profile (either by creation of a distinct 
bulge component, or simply by steepening the light profile) 
but manages to retain a significant fraction of the system's 
original angular momentum in preserving a disk component. 
In the context of galaxy merging, such systems are a relatively 
natural outcome of merging between disk galaxies with even 
modest gas fractions, where higher mass ratio minor mergers 
(or mergers between more gas-rich systems) lead to progres- 
sively more disk-dominated remnants (see, e.g., iNaab et aTl 
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Fig. 1 1. — F160W Postage stamps of quiescent galaxies witli u < 2 and 
stellar surface densities between lO' and IO'^MqIcpc"^. These galaxies 
should be star-forming based on their Sersic indices, but they are in fact qui- 
escent. The postage stamps, within each class, are ordered by stellar mass 
(ordered left to right). At all redshifts, all postage stamps are 40 physical kpc 
on a side, and are scaled to a constant 'stellar mass density' (total intensity is 
scaled to total stellar mass, meaning that if the stellar M/L is constant over the 
face of the galaxy this postage stamp should reflect the stellar mass density) 
and are displayed using asin h scaling (lin ear at low intensity, and logarithmic 
at higher surface brightness: ILupton et alJ1999l) . 



l2(MlHopkins et al.ll2()()9llHoffman et alJIMoh . 

Secondly, a small fraction of quiescent systems a ppear 
to have little in the way of a bulge comp onent dBell 



20081: IStockton et alj 20Q4; McGrath et al. 2008; Bundv et al 



20101); the examples in McGr ath et al.i (i2008i) have been very 
carefully documented. At least a few of the examples in Fig. 



[TT] appear to have genuinely low-concentration light profiles, 
and at least one is an inclined disk. Given that such systems 
in the local universe appear to all be satellites in high mass 
(group or cluster) halos (Bell 2008), it will be interesting to 
explore the environments of such galaxies, as a function of 
redshift, and help to elucidate the extent to which disk-only 
quiescent galaxies are the products of stripping of their cold 
gas content by hot gas in a deep potential well (in this context, 
it is worth noting that the z 1 .5 systems in iMcGrath et al.l 
12008 were chosen to be in fields near radio-loud z ~ 1.5 
QSOs, and may in fact reside in overdensities). In this con- 
text, we note that R. Bassett et al. (in preparation) have ex- 
plored this issue in a narrow redshift slice with the dataset 
used in this paper, focusing on an overdensity at z ^ 1.62, 
finding weak evidence of an increase in the fraction of qui- 
escent galaxies with low Sersic indices near the largest over- 
density in the UDS. Their result, although it must be explored 
with larger samples, appears consistent with the notion that 
disk-only quiescent galaxies should be best interpreted as be- 
ing gas-free owing to external (e.g., environmental) mecha- 
nisms, not internal processes. 

5.3. On the nature of high Sersic index star-forming galaxies 

Figs. [8] and [TO] show a population of star-forming galax- 
ies with high Sersic index. We show examples of such sys- 
tems in Fig. [12] in 5 different redshift bins. One can see that 
these high Sersic index, star-forming systems fall broadly into 
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n>3 star forming galaxies 



• • • • 

1.3<z<1.5 

n>3 star forming galaxies 




0.9<z<1.3 

n>3 star forming galaxies 

0.6<z<0.9 » 



Fig. 12. — F160W Postage stamps of star-forming galaxies with n> 3 and 
stellar surface densities between lO' and IO'^Mq kpc"^. The postage stamps, 
within each class, are ordered by stellar mass (ordered left to right). At all 
redshifts, all postage stamps are 40 physical kpc on a side, and are scaled 
to a constant 'stellar mass density' (total intensity is scaled to total stellar 
mass, meaning that if the stellar M/L is constant over the face of the galaxy 
this postage stamp should reflect the stellar mass density) and are displayed 
using asin h scaling (linear at l ow intensity, and logarithmic at higher surface 
brightness JLupton et al.ll999D . 

two types of systems. At lower redshifts, there are some 
clear examples of disks with very prominent bulges; such 
galaxies appear to be the star-forming counterparts of the 
quiescent disks discussed above, and appear to be systems 
that have either retained, or 're-grown', a substantial disk 
during/after the bulge formation process jBa ugh et al. 1996; 
iKannappa n et alj 120091; iBenson & Devereux 2010). At all 
redshifts, there are what appear to be genuinely spheroidal 
galaxies, but in many cases with significant asymmetries, and 
in some cases evidence of interactions (e.g., nearby compan- 
ions, tidal tails). It will be interesting to examine such sys- 
tems in the future to understand the degree of over lap be- 
tween these systems and 'blue sphero ids' (Mena nteau et al. 
[2001; HauBler 2007; IKannappan et al.1 2009; G vorv & Bell 
|2010), post-starburst galaxies ("Vergani et al. 201 or or ongo- 
ing 'bulge-building' starbursts (Wuyts et al. 201 lb) and to ask 
if the properties of these galaxies are more consistent with an 
interpretation as a quiescent galaxy in formation, a galaxy in 
the midst of bulge building through instabilities in the stel- 
lar disk and cold inflowin g gas (Dek el, Sari & Ceverino„2009; 
ICeverino. Dekel & Bourn aud 2010 ), or a galaxy in the early 
stages of re-growing a disk (iKannappan et alJ[2009^ ■ 

5.4. Future steps 

There are a number of possible future directions for pur- 
suing this line of study further A clear improvement path 
involves the use of more accurate photometric spectroscopic 
redshifts, which helps improve both the accuracy and relia- 



bility of redshift and stellar population estimates, and opens 
up possibilities to measure low-resolution spectral parame- 
ters like the 4000A break sti'ength (e.g.. Whit aker et aUlIoTll; 
l&iek et al. 201 1; van Dokkum et al. 201 1). Further sharpen- 
ing the analysis would be possible with the addition of reliable 
bulge/disk decompositions: the hypothesis that the existence 
of a prominent bulge is necessary but not sufficient to shut off 
star formation can be better tested, and additional parameters 
can be explored for possible relevance (e.g., bulge/total ratio, 
bulge Sersic index, or bulge mass; see, e.g.. iDrorv & Fishen 
1200 71 for an exploration of some of these issues with a low- 
redshift sample, or V. Bruce et al. in preparation). More 
detailed exploration of the outliers at all redshifts (the qui- 
escent lower Sersic index systems, or star-forming systems 
with high Sersic index) may help to illuminate the processes 
that remove (or keep out) cold gas from galaxies. Finally, 
increasing the number statistics at z < 1 (already underway 
with e.g., the GEMS, AEGIS or CO SMOS datasets; lRix et alj 
l2004llDavis et allllOOTllScoville et a l. 2007) andatz> Iwifli 
the fuU five- field coverage of CANDELS (iGrogin et al.ll2011l; 
iKoekemoer et al. 201 1) is likely to prove useful. 

6. CONCLUSIONS 

Motivated by a desire to empirically explore the evolution- 
ary factors that lead to a lack of star formation in galaxies, 
we have explored the structures and star formation activity 
of the galaxy population using the HST/WFC3 F160W imag- 
ing data from CANDELS in the UDS field. We used pub- 
lic pho tometry and photometric redshifts from lWilliams et aT] 
1*2009), and determined rest-frame absolute magnitudes and 
stellar masses using our own stellar population model fits. We 
supplement this with public 24/im data from SpUDS, and sep- 
arate galaxies into quiescent and star-forming using a combi- 
nation of optical-near-infrared colors and 24/im information. 
For structural information, we use parametric fits of a single 
Sersic profile to the WFC3 F160W imaging data. We then 
combine these data for z > 0.6 galaxies with a sample with 
similar parameters from the SDSS, to create a z < 0.05 com- 
parison sample. We then proceed to explore the evolution of 
the > 3 X 1O'°M0 galaxy population from z = 2.2 to the 
present day. 

We first visualize the evolution of the galaxy population 
over the last 10 Gyr by normalizing the sample to a 'fixed' 
comoving volume of 10^ Mpc^. In agreement with a large 
number of other works, we find that the number density of 
galaxies with stellar mass in excess of 3 x 1O'"M0 increases 
approximately five-fold from z ^ 2 to the present day, and that 
the number density of quiescent galaxies increases yet more 
rapidly. Furthermore, examining the properties of the quies- 
cent galaxy population, we find that the vast majority of those 
quiescent galaxies have high Sersic indices (and, at z > 1, high 
inferred velocity dispersion and surface density), a sign that 
they have a prominent bulge component. The growth of the 
quiescent galaxy population appears to be intimately linked 
with the growth of galaxies with prominent bulges. 

Given the rapid evolution of the quiescent galaxy popula- 
tion, we proceed to explore the strength of the correlation be- 
tween quiescence and four galaxy/structural parameters: stel- 
lar mass, 'velocity dispersion' M/R oc a'^ with a Sersic index- 
dependent scaling, stellar surface density and Sersic index. At 
all redshifts z < 2.2, stellar mass correlates poorly with qui- 
escence. It is possible, bearing in mind the inferred < 0.2 dex 
scatter between stellar mass and dark matter halo mass at z ~ 
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dMore et alj|2009l) . that the weakness of this correlation indi- 
cates that halo mass alone is not the ma in determinant of qui- 
escence (but see also lWake et al.l2012al and lWake et al.ll2012bl 
for arguments that velocity dispersion correlates with halo 
mass better than stellar mass, at least at high stellar mass). 

At z < 0.05, we find that quiescence correlates much more 
strongly with Sersic index than either velocity dispersion or 
surface density. At z > 0.6, we find that velocity dispersion, 
surface density and Sersic index correlate well with quies- 
cence, where the correlation of Sersic index with quiescence 
appears marginally stronger 

All correlations have substantial scatter, however Many 
quiescent systems have prominent disks, although the vast 
majority of quiescent galaxies with disks also have prominent 
bulges. A very small fraction of quiescent galaxies appear to 
be bulgeless disks. In the local Universe they are all satel- 
lite galaxies in galaxy groups/clusters; at higher redshifts we 
did not explore environmental variables for lack of dynamic 
range in galaxy environments (see R. Bassett et al. in prepa- 
ration for first steps in this direction using the CANDELS 
dataset). Star-forming systems with high n, M, S and a' are 
not particularly uncommon; at high redshifts they appear to 
be genuinely compact with high n, and often show asymme- 
tries or signatures of tidal interactions (one may wish to asso- 
ciate these with the possible remnants of gas-rich dissipational 
galaxy interactions/mergers), and at lower redshifts there is a 
mix of similar systems and composite bulge-disk star-forming 
systems (with large bulges). 

At z < 0.05, Bell (2008) concluded that a prominent bulge 
(and by association, a supermassive black hole) was a neces- 
sary but not sufficient condition for a galaxy to turn off its own 
star formation on galaxy-wide scales (all quiescent galaxies 
in the centers of their own halos had prominent bulges, but 
not all galaxies with bulges lack star formation). This obser- 



vational association is qualitatively consistent with the AGN 
feedback paradigm (no supermassive black hole, no ability to 
shut off star formation). While there is clearly scope for fur- 
ther investigation of the drivers of quiescence, the evidence 
assembled here appears to be consistent with this proposition 
to z < 2.2, a time interval of more than lOGyr. 
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APPENDIX 

APPENDIX A: NORMALIZING THE SAMPLE TO A CONSTANT COMOVING VOLUME 

In Fig.m we Monte Carlo subsampled the galaxy population of the UDS and the SDSS to an equivalent comoving volume of 
lO^Mpc^. There are two ways to do this; simply rescaling the sample by 10^Mpc^/y(z), where V(z) is the comoving volume 



20 



Bell et al. 




+ <> 



I 200 - 





^ 100- 



o 



CD 



50 - 



E 







0.0 0.5 1.0 1.5 2.0 



z 



Fig . 13. — The number of galaxies with stellar masses in exce ss of 3 X lO^^M(r, in a v olume of lO^Mpc^ measured from the UDS (diamonds, with the dominant 
sample variance uncertainties as estimated using the method of lSomerville et al."2004': see also lMoster e t al.''201 V). and the number of galaxies with masses in 
excess of 3 X 1()"'Mq (crosses) infeired from the stellar mass functions offllbert et al. (2010, z > 0.2) and Bell et al. (2003, local values). The number of galaxies 
with Mf > 3 X IO'^Mq in the UDS is consistent to within the significant sample variance and systematic stellar mass uncertainties with those of larger fields. We 
choose in this paper to adjust the number of galaxies for display in Fig.|4]to fit the straight line in this figure, to take out the first-order effects of sample variance 
on the results (while noting that doing so makes no important difference to our conclusions; the evolution of the population is not subtle). 

of that redshift bin (leaving one susceptible to the first-order effects of sample variance), or by rescaling the sample in a given 
redshift bin to have the number of galaxies expected at that redshift above that mass limit using mass functions derived from 
much larger surveys (canceling out number density variation from sample variance but leaving behind any systematic variation 
in galaxy properties that are a function of the variation in the average environment in this light cone as a function of redshift). 

The two approaches are compared in Fig. [13] Diamonds with error bars denote the number of galaxies with stellar masses 
m excess of 3 X 10'°Mq in the UD S, with error bars denoting the expected degree of sample variance in that bin following the 
method of 'Som erville et afl (120041) . Crosses at z > 0.2 show the number of galaxies with stellar mass in excess of 3 x 1O'*'M0 
expected from the stellar mass functions of Ilbert et al. (2010). The cross at low redshift shows the number of galaxies with 
>3x IO'^Mq from the mass function of Bell et al. (2003). The number of galaxies observed in the UDS is broadly consistent 
with, or perhaps somewhat larger than, the number of galaxies expected from larger surveys, given the substantial sample variance 
uncertainties. There are systematic uncertainties also on the crosses; choosing different star formation histories for constructing 
stellar masses, stellar population models, etc., can give more than a factor of two variation in stellar masses that translates into 
around 50% in number density uncertainty. For Fig. |4] we chose to rescale the number of galaxies to the smoothly-varying 
number of galaxies given by the line shown in Fig.[T2](approximately corresponding to scaling the number of galaxies to larger 
cosmological surveys). 

None of the results shown in Fig. |4] depend in any important way on the choice of this scaling method; the evolution of the 
galaxy population from z ~ 2 to the present day is not subtle and is robust to even the significant systematic uncertainties inherent 
to mass function analyses. While this scaling in the number of galaxies largely counteracts the worst effects of sample variance, 
if the properties of galaxies depend strongly on environment there will be a second-order difference between the properties of 
galaxies sampled from a globally-underdense vs. a globally-overdense volum e. In the absence of en ough environmental dynamic 
range to robustly measure its effect within CANDELS at the current time (see lPapovich et all2012l and R. Bassett et al., in prep., 
for a discussion of SFR, galaxy size and structure trends as a function of environment using a known z= 1 .62 galaxy cluster in 
the CANDELS UDS coverage), we made no attempt to correct for this second-order effect in this paper. 



